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Abstract 

We propose a numerical method to evaluate the upper critical dimension dc of random percolation 
clusters in Erdos-Renyi networks and in scale-free networks with degree distribution V{k) ~ k^'^, 
where k is the degree of a node and A is the broadness of the degree distribution. Our results 
report the theoretical prediction, dc = 2(A — 1)/(A — 3) for scale-free networks with 3 < A < 4 and 
dc = 6 for Erdos-Renyi networks and scale- free networks with A > 4. When the removal of nodes 
is not random but targeted on removing the highest degree nodes we obtain = 6 for all A > 2. 
Our method also yields a better numerical evaluation of the critical percolation threshold, Pc, for 
scale-free networks. Our results suggest that the finite size effects increases when A approaches 3 
from above.. 

PACS numbers: 
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Recently much attention has been focused on the topic of complex networks, which char- 
acterize many natural and man-made systems, such as the Internet, airline transport system, 
power grid infrastructures, and the world wide web (WWW) 

HQS 

, |j| . Many studies on 

these systems reveal a common power law degree distribution, V{k) ~ k^^ with k > /Cmim 
where k is the degree of a node, A is the exponent quantifying the broadness of the de- 
gree distribution 5| and kmm is the minimum degree. Networks with power law degree 
distribution are called scale-free (SF) networks. The power law degree distribution repre- 
sents topological heterogeneity of the degree in SF networks resulting in the existence of 
hubs that connect significant fraction of nodes. In this sense, the well studied Erdos-Renyi 
(ER) networks are homogeneous and can be represented by a characteristic degree 

(k), the average degree of a node, while SF networks are heterogeneous and do not have a 
characteristic degree. 

The embedded dimension of ER and SF networks can be regarded as infinite {d = oo) 
since the number of nodes within a given "distance" increases exponentially with the dis- 
tance compared to an Euclidean d dimensional lattice network where the number of nodes 
within a distance L scales as L'^. Percolation theory is a powerful tool to describe a large 
number of systems in nature such as porous and amorphous materials, random resistor net- 
works, polymerization process and epidemic spreading and immunization in networks [o], [l^ . 
Percolation theory study the topology of a network of nodes resulting from removal of a 
fraction q = 1 — p of nodes (or links) from the system. It is found that in general there exists 
a critical phase transition at p = pc, where pc is the critical percolation threshold. Above pc, 
most of the nodes (order A^) are connected, while below pc the network collapses into small 
clusters of sizes of order InA^. For lattices in c? > 6, the nodes, in the percolation cluster, 
do not have spatial constraints and therefore all percolation exponents remain the same and 
the system behavior can be described by mean field theory [o], [lo|. This is because &t dc = 6 
the spatial constraints on the percolation clusters become irrelevant and each shortest path 
between two nodes in the percolation cluster at criticality can be considered as a random 
walk. The critical dimension dc above which the critical exponents of percolation become 
the same as in mean field theory is called the upper critical dimension (UCD). It is well 
known that the UCD for percolation in d-dimensional lattices is 6. Studies of percolation in 
ER networks, yield the same critical exponents as in mean-field values of regular percolation 
in infinite dimensions. This is because in ER networks spatial constraints do not appear and 
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the symmetry is almost the same as in Euchdean lattices, i.e., there is a typical number of 
links per node. However, SF networks with 2 < A < 4 have different critical exponents than 
ER networks 
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121 ]. The regular mean-field exponents are recovered only for SF networks 
with A > 4. This is due to the fact that for the classical mean field one needs two conditions 
(a) no spatial constraint (b) translational symmetry, meaning that all nodes have similar 
neighborhood. The second condition does not apply for SF networks with A < 4 due to 



the broad degree distribution and thus we expect a new type of mean field exponents {4]. 
Indeed, for SF networks with 3 < A < 4, the UCD was shown to be [l^ : 

Thus, dc is larger than 6 and for A — >• 3, dr — >■ 00. When scale- free networks are embedded 



in a regular Euclidean lattice [13|, ll5ll) the value of dc tells us above which dimension 
the percolation clusters will not be affected by the spatial constraints and therefore the 
percolation exponents will be the same as for infinite dimension. Thus, it is reasonable that 
when A is smaller, the network is more complex (due to bigger hubs) and a higher upper 
critical dimension is expected. However, Eq. ([1]), that was shown analytically to be valid 
for N ^ 00 was never verified or tested numerically. It is also interesting to determine the 
range of values where the results of Eq. ([T]) can be observed. Here we propose a numerical 
method to measure the value of dc for ER and SF networks with A > 3 |16|. 

Finite-size scaling arguments in d-dimensional lattice networks predict [ol, [l^ that the 
critical threshold Pc{L) approaches Pc = Pc(oo) via, 

p,(L)-p,(oo)~L-V^ (2) 

where L is the linear lattice size and u is the correlation critical exponent. Eq. ([2]) for lattices 
can be generalized to networks of nodes via the relation L'^ = N, i.e., Pc{N) — pc{oo) ~ 
Since networks can be regarded as embedded in infinite dimension and since above 
dc all exponents are the same, we replace d by dc, 

Pc{N) - pc{oo) ~ N-^'"^'' = iV^®. (3) 

For ER and SF networks with A > 4, we have dc = Q and u = 1/2, thus from Eq. ([3]) follows, 

p,(iV)-p,(oo)~iV-V3. (4) 
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For SF networks with 3 < A < 4, we have z/ = 1/2 and substituting Eq. ([T]) in Eq. ([3]), it 
yield, 

p,(iV)-p,(oo)~iV(3-^)/(^-i). (5) 

In this paper we use Eq. to measure = 2/dc from which we can evaluate dc- To 
measure 9, using the finite size scaling of Eq. ([3]), we have to compute the dependence of 
the percolation threshold, pc{N), of ER and SF networks on the system size A^. To calculate 
Pc{N), we apply the second largest cluster method [o, Q, which is based on determining 
Pc{N) by measuring the value of pc at the maximum value of the average size of the second 
largest cluster, (S'2). It is known that {S2) has a sharp peak as a function of p at Pc 0, 0. To 
detect this peak we perform a Gaussian fit around the peak and estimate the peak position 
which is pc{N) [3]. 

To improve the speed of the simulations, we implement the fast Monte Carlo algorithm 
for percolation proposed by Newman and Ziff [l^ . Basically, for each realization, we prepare 
one instance of N nodes network with the desired structure as the reference network. Then 
we prepare another set of nodes with no links as our target network. Because we want 
to know the size of the 2nd largest cluster instead of the largest one, we use a list which 
keeps track of all the clusters in descending order according to their sizes, which in the 
beginning is a list of clusters of size one. As we choose the links in random order from the 
reference network and make the connection in the target network, we update the list of the 
cluster size but always keep them in descending order. The concentration value, p, of each 
newly connected link is calculated by the number of links after adding this link in the target 
network divided by the total number of links in the reference network. We record S'2 in the 
following way. First, we make 1000 bins between and 1. When each link is connected, we 
record S'2 at the concentration value p of this newly connected link. After many realizations, 
we take the average of S2 for each bin. 

Figure [Dj^a) shows (S2) as a function of p, for two different system sizes of ER networks 
with (k) = 4. The position of the peak, obtained by fitting the peak with a Gaussian 
function, yields Pc{N). Figure [I](b) shows Pc{N) as a function A^. Using Pc(oo) = l/{k) = 
0.25 0, Q, the fitting of Eq. ([3]) gives the exponent 6 = 0.328 ± 0.003, very close to the 
theoretical prediction for ER, = 1/3, Eq. (jl]). We performed the same simulations for ER 
with other average degrees, (k) = 5 and 6, and obtained similar results for 6. 
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To determine Pc{oo) for random SF networks, we use the exact analytical results |19| . 

Pc(oo) = — (6) 

Here kq = (^o)/(^o) is computed from the original degree distribution {V{kQ)) for which the 
network is constructed. However, the way to compute the value of kq is strongly affected by 
the algorithm of generating the SF network as explained below. 

To generate SF networks with power law exponent A, we use the MoUoy-Reed algo- 



rithm 
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2l\ . We first generate a series of random real number u satisfying the distribution 



V{u) = cu~^, where c = (A — l)//i;^7n is the normalization factor. Next we truncate the real 
number u to be an integer number k, which we assume to be the degree of a node. We make 
k copies of each node according to its degree and randomly choose two nodes and connect 
them by a link. Notice that the process of truncating the real number u to be an integer 
number k which is the degree of a node actually slightly changes the degree distribution 
because any real number n < u < n + 1, where n is an integer number, will be truncated to 
be equal n. Thus, the actual degree distribution we obtain using this algorithm is 



fc+i ^ ^ 1 
k 



V{k) = I cu-^du = TTZX (^'"^ + 1)'"^). (7) 



mm 



We use Eq. ([7]) to compute and Pc(oo) defined in Eq. IQ. Table [T] shows the calculated 
results of Pc(oo) for several values of A. 

We calculate (5*2) for SF networks for different values of A and and compute Pc{N) by 
fitting with a Gaussian function near the peak of (5*2) as for ER networks. Using the values 
of Pc{oo) for SF networks displayed in Table [H we obtain B by a power law fitting with 
Eq. ([3]) as shown in Fig. [21 As we can see for A = 4.5, 3.85 and 3.75 we obtain quite good 
agreement with the theoretical values. However for A = 3.65 and 3.5, the values of G become 
better when fitting only the last several points (largest A^) and still have large deviations 
from their theoretical values. This strong finite size effect isprobably since for A — > 3 the 
largest percolation cluster at the criticality becomes smaller 22] . Thus, we expect that as A^ 
increase, the exponent Q{N) obtained by simulations should approach the theoretical value 
of G of Eq. ([5]). To better estimate G we assume finite size corrections to scaling for Eq. (Q, 
i.e., 

p,{N) - p,(oo) ~ Ar-e(l + N~-). (8) 
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Thus, the actual 0(A^) obtained from simulation is the succesive slopes, 



e(iV) = -d\n{p,{N) - p,ioo))/di\nN), (9) 

from which we can see that B(A^) approaches B as a power law, 

Q{N) - ~ A^-^. (10) 

Indeed, Fig. [3] shows the exponent G(A^) as a function of A^~^ for A = 3.5 and 3.65. Fig- 
ure[3]^a) shows that for A = 3.5 and x = 0.11, we obtain a straight line and Q{N) approaches 
0.2 as — > oo, consistent with the theoretical value of O (Table [11). Fig. [3](b) shows, for 
A = 3.65 and x = 0.13, O(A^) is again a straight line that approaches 0.245 for —>■ oo, 
consistent with the theory. 

Next we estimate the value of dc for SF network under targeted attack on the largest 



degree nodes 
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251]. For this case since the hubs are removed we expect that for all 
\ > 2, dc will be the same as for ER, i.e., dc = 6. In Fig. HJ we plot Pc{N) — Pc{oo) for SF 
with A = 2.5 under targeted attack. Indeed from Eq. ([3]) by changing Pc{oo) and fitting the 
best straight line in log-log plot, we obtain 9 ~ 0.33, i.e., dc ~ 6, as expected. 

Further supports of the analytical aproach, we evaluat by simulations V{s), the prob- 
ability distribution of the cluster sizes at pc{N), which should follow a power law for SF 



networks 



V{s) ~ s-^ = 2 < A < 4. (11) 



Figure [5] shows the simulations results for SF networks A = 3.5. The dashed line is the 
reference line with slope —2.67, which is the theoretical value of r from Eq. (I5l), showing 
good agreement between theory and simulations. 
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FIG. 1: (a) The average size of the 2nd largest cluster, (5*2), as a function of the concentration, 
p, of links present in the ER networks. The typical number of realizations for each curve is 10^. 
(b) Log-log plot pc{N) — pc{oo) as a function of N, where Pc{oo) = l/{k) = 0.25 for ER with 
{k) = 4. 

[26] L. K. Gallos, et al., Phys. Rev. Lett. 94, 188701 (2000). 



A 


Pc(oo) 


Theoretical 


Numerical @ 


3.50 


0.2039 


0.200 


0.234 


3.65 


0.2574 


0.245 


0.260 


3.75 


0.2911 


0.273 


0.275 


3.85 


0.3234 


0.298 


0.284 


4.50 


0.5009 


1/3 


0.326 


ER {{k) = 4) 


0.25 


1/3 


0.328 



TABLE 1: The main results for SF and ER networks. The critical percolation threshold pdoo) 
indicates the numerical value calculated according to Eqs. ([6]) and ([7]). Theoretical Q is the 
theoretical prediction of Q (from Eqs. ([T])) and ^ and numerical is the numerical value we 
obtained from simulations. The SF networks were generated with fcmin = 2. 
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(e) X=4.5 
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(d) X=3.85 
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FIG. 2: Log-log plots of Pc{N) — Pc(oo) as a function of N for SF networks with /cmin = 2 and 
different value of A. The dashed line is the reference line with indicated slope. 
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FIG. 3: The exponent O(A^) as a function of N^'-^ for SF networks with k^nhi = 2 and different 
value of A: (a) A = 3.5, where x ^ 0.11; and (b) A = 3.65, where x ^ 0.13. The theoretical values 
9(oo) = 0.2 (A = 3.5) and 9(cxd) = 0.245 (A = 3.65), are consistent with the asymptotic values of 
9 obtained for A?^ — oo. 




FIG. 4: Log-log plot of Pc{N) ~ Pc{oo) as a function of N for SF networks with A = 2.5, /cmin = 2 
for a targeted attack. The dashed line is the best fit with slope —0.33. Since we do not have a good 
estimation for pc{oo), we modified Pc{oo) to get the best straight line in log-log plot, Pc{oo) = 0.23 
(o), Pc{oo) = 0.25 (□) and Pc{oo) = 0.26 (O). When Pc{N) — Pc{oo) is linear (dashed line) in the 
log-log plot, the slope yields the exponent 9 !=a 0.33 i.e., dc = 6. 
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FIG. 5: The probability distribution of the cluster sizes at Pc{N) for TV = 2048 (o) and N = 16384 
(□). The dashed line is the reference line with slope —2.67. 
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